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I examine electron-phonon mediated superconductivity in the intermediate coupling and phonon 
frequency regime of the quasi-2D Holstein model. I use an extended Migdal-Eliashberg theory which 
includes vertex corrections and spatial fluctuations. I find a d-wave superconducting state that is 
unique close to half-filling. The order parameter undergoes a transition to s-wave superconductivity 
on increasing filling. I explain how the inclusion of both vertex corrections and spatial fluctuations is 
essential for the prediction of a d-wave order parameter. I then discuss the effects of a large Coulomb 
pseudopotential on the superconductivity (such as is found in contemporary superconducting ma- 
terials like the cuprates), which results in the destruction of the s-wave states, while leaving the 
d-wave states unmodified. Published as: Phys. Rev. B 73, 060503(R) (2006) 
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The discovery of high transition temperatures and a 
d-wave order parameter in the cuprate superconductors 
are remarkable results and have serious implications for 
the theory of superconductivity. The presence of large 
Coulomb interactions in the cuprates which have the po- 
tential to destroy conventional s-wave BCS states has 
prompted the search for new mechanisms that can give 
rise to superconductivity. However, electron-phonon me- 
diated superconductivity is still not well understood, es- 
pecialy in lower dimensional systems. In particular, the 
electron-phonon problem is particularly difficult at inter- 
mediate couplings with large phonon frequency (such as 
found in the cuprates) and the electron-phonon mech- 
anism cannot be fully ruled out. It is therefore of 
paramount importance to develop new theories to under- 
stand electron-phonon mediated superconductivity away 
from the BCS limit. 

The assumption that electron-phonon interactions can- 
not lead to high transition temperatures and unusual or- 
der parameters was made on the basis of calculations 
from BCS theory, which is a very-weak-coupling mean- 
field theory (although of course highly successful for 
pre-1980s superconductors) 1 . In the presence of strong 
Coulomb interaction, the BCS s-wave transition temper- 
ature is vastly reduced. However, the recent measure- 
ment of large couplings between electrons and the lat- 
tice in the cuprate superconductors means that exten- 
sions to the conventional theories of superconductivity 
are required2i^. In particular, low dimensionality, in- 
termediate dimensionless coupling constants of ~ 1 and 
large and active phonon frequencies of ~ 75meV mean 
that BCS or the more advanced Migdal-Eliashberg (ME) 
theory cannot be applied. In fact, the large coupling con- 
stant and a propensity for strong renormalization in 2D 
systems, indicate that the bare unrenormalized phonon 
frequency could be several times greater than the mea- 
sured 75 meV— . 

Here I apply the dynamical cluster approxima- 
tion (DCA) to introduce a fully self-consistent 
momentum-dependent self-energy to the electron-phonon 



problem^i^. Short ranged spatial fluctuations and low- 
est order vertex corrections are included, allowing the 
sequence of phonon absorption and emission to be re- 
ordered once. In particular, the theory used here is 
second order in the effective electron-electron coupling 
U = —g 2 /AluiQ, which provides the correct weak coupling 
limit from small to large phonon frequencies^. In this 
paper, I include symmetry broken states in the anoma- 
lous self energy to investigate unconventional order pa- 
rameters such as d-wave. No assumptions are made in 
advance about the form of the order parameter. 

DCA 6,8,9 is an extension to the dynamical mean-field 
theory for the study of low dimensional systems. To ap- 
ply the DCA, the Brillouin zone is divided into Nc sub- 
zones within which the self-energy is assumed to be mo- 
mentum independent, and cluster Green functions are 
determined by averaging over the momentum states in 
each subzone. This leads to spatial fluctuations with 
characteristic range, Nc /D . In this paper, N c — 4 is used 
throughout. This puts an upper bound on the strength of 
the superconductivity, which is expected to be reduced 
in larger cluster sizes^S. To examine superconducting 
states, DCA is extended within the Nambu formalism 7-8 . 
Green functions and self-energies are described by 2 x 2 
matrices, with off diagonal terms relating to the super- 
conducting states. The self-consistent condition is: 



G(K,iuj n ) = / rfe 



A(e)(C(K. t ^ ra )-6) 
|C(K;, iuj n ) - e\ 2 + 4>(Ki,iuj n ) 2 



(1) 



F(K, iu n ) 



de 



"Di(e)4>(Ki,iu} n ) 



|C(K i ,^„)-e| 2 + 0(K J ,i W „) 2 

(2) 

where C(K.j, iuj n ) — iuj n + fi — £(Kj, iu n ), fi is the chem- 
ical potential, tu n are the Fermionic Matsubara frequen- 
cies, </>(K, iuj) is the anomalous self energy and £(K, iu>) 
is the normal self energy. G(K, iuj n ) must obey the lat- 
tice symmetry. In contrast, it is only |.F(K, iw n )\ which 
is constrained by this condition, since <p is squared in 
the denominator of Eqn. [TJ Therefore the sign of <f> 
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FIG. 1: Diagrammatic representation of the current approx- 
imation. Series (a) represents the vertex-neglected theory 
which corresponds to the Migdal-Eliashberg approach, valid 
when the phonon energy u>o and electron-phonon coupling 
U are small compared to the Fermi energy. Series (b) repre- 
sents additional diagrams for the vertex corrected theory. The 
phonon self energies are labeled with II, and E denotes the 
electron self-energies. Lines represent the full electron Green 
function and wavy lines the full phonon Green function. 



can change. For instance, if the anomalous self energy 
has the rotational symmetry (f)(w, 0) = — <^>(0, 7r), the 
on-diagonal Green function, which represents the elec- 
tron propagation retains the correct lattice symmetry 
G(w, 0) = G(0,7r). Therefore, only inversion symmetry 
is required of the anomalous Green function representing 
superconducting pairs and the anomalous self energy. 

Here I examine the Holstein model 11 of electron- 
phonon interactions. It treats phonons as nuclei vibrat- 
ing in a time-averaged harmonic potential (representing 
the interactions between all nuclei), i.e. only one fre- 
quency luq is considered. The phonons couple to the local 
electron density via a momentum-independent coupling 
constant g~. 
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The first term in this Hamiltonian represents hopping of 
electrons between neighboring sites and has a dispersion 
£k = — 2t Y^h=i cos(Jbj). The second term couples the lo- 
cal ion displacement, to the local electron density. The 
last term is the bare phonon Hamiltonian, i.e. a sim- 
ple harmonic oscillator. The creation and annihilation of 
electrons is represented by c\(ci), pi is the ion momen- 
tum and M the ion mass. The effective electron-electron 
interaction is, 
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where, lo s = 2irsT, s is an integer and U = —g 2 /Mu>Q 
represents the magnitude of the effective electron- 
electron coupling. D = 2 with t = 0.25, resulting in 
a non-interacting band width W = 2. A small intcrpla- 
nar hopping t± = 0.01 is included. This is necessary to 
stabilise superconductivity, which is not permitted in a 
pure 2D system^. 



Perturbation theory in the effective electron-electron 
interaction (Fig. [T|) is applied to second order in U, us- 
ing a skeleton expansion. The electron self-energy has 
two terms, Eme(w,K) neglects vertex corrections (Fig. 
QJa)), an d £vc(w, K) corresponds to the vertex corrected 
case (Fig. [1Kb)). Hme(w,K) and Hvc(^,K) correspond 
to the equivalent phonon self energies. At large phonon 
frequencies, all second order diagrams including Y,yc are 
essential for the correct description of the weak coupling 
limit. 

The phonon propagator D(z, K) is calculated from, 
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and the Green function from equations [T] and [2 S = 
S ME + S vc and n = + n-vc • Details of the transla- 
tion of the diagrams in Fig. [T]and the iteration procedure 
can be found in Ref. 0- Calculations are carried out along 
the Matsubara axis, with sufficient Matsubara points for 
an accurate calculation. The equations were iterated un- 
til the normal and anomalous sclf-cncrgics converged to 
an accuracy of approximately 1 part in 10 3 . 

Since the anomalous Green function is proportional 
to the anomalous self energy, initializing the problem 
with the non-interacting Green function leads to a non- 
superconducting (normal) state. A constant supercon- 
ducting field with d-wave symmetry was applied to the 
system to induce superconductivity. The external field 
was then completely removed. Iteration continued with- 
out the field until convergence. This solution was then 
used to initialize self-consistency for other similar val- 
ues of the parameters. The symmetry conditions used in 
Refsiand0 have been relaxed to reflect the additional 
breaking of the anomalous lattice symmetry in the <i-wave 
state. This does not affect the normal state Green func- 
tion, but does affect the anomalous state Green function. 

In Fig. [21 the anomalous self energy is examined for 
n = 1.0 (half- filling) . The striking feature is that sta- 
ble d-wave superconductivity is found. This is mani- 
fested through a change in sign of the anomalous self 
energy, which is negative at the (it, 0) point and positive 
at the (0, 7r) point. The electron Green function (equa- 
tion [lj depends on (f> 2 , so causality and lattice symmetry 
are maintained. Since the gap function <p(iuj n ) / Z (iu) n ) 
is directly proportional to (j)(iu> n ), and Z(iu) n , K^o)) = 
Z(iu} n , K(o )7r )), then the sign of the order parameter i.e. 
the sign of the superconducting gap changes under 90° 
rotation. Z(iu n ) = 1 — T,(iuj n ) j 'iui n . 

Figure fS] shows the variation of superconducting pair- 
ing across the Brillouin zone. rt s (k) = T£ n F(iw„,k). 
U = 0.6, ujq = 0.4,71 = 1 and T = 0.005. The d-wave 
order can be seen very clearly. The largest anomalous 
densities are at the (71", 0) and (0, 7t) points, with a node 
situated at the (ir/2, n/2) point and a sign change on 90° 
rotation. Pairing clearly occurs between electrons close 
to the Fermi surface. 

So far, the model has been analyzed at half filling. 
Figure 2] demonstrates the evolution of the order param- 
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FIG. 2: Anomalous self-energy at half-filling. The anomalous 
self energy is real. It is clear that 0(tt, 0) = — 0(0, 7r). This 
is characteristic of d-wave order. Similarly, the electron self 
energy has the correct lattice symmetry E(-7r,0) = B(0, 7r), 
which was not imposed from the outset. The gap function is 
related to the anomalous self energy via cj>(iui n )/Z(iu} n ). 
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FIG. 3: Variation of superconducting (anomalous) pairing 
density across the Brillouin zone. n s (k) ~ T^2 F(iui n ,'k). 
U = 0.6, lo = 0.4, n = 1 and T = 0.005. The d"wave order 
can be seen very clearly, with a change in sign on 90° rota- 
tion and a node situated at the (n/2, n/2) point. The largest 
anomalous (superconducting) densities are at the (it, 0) and 
(0, 7r) points. 



eter as the number of holes is first increased, and then 
decreased. The total magnitude of the anomalous den- 
sity, n s — J2k l n s(Kj)| is examined. When the number 
of holes is increased, stable d-wave order persists to a 
filling of n — 1.18, while decreasing monotonically. At 
the critical point, there is a spontaneous transition to s- 
wave order. Starting from a high filling, and reducing the 
number of holes, there is a spontaneous transition from 
s to d-wave order at n = 1.04. There is therefore hys- 
teresis associated with the self-consistent solution. It is 
reassuring that the d-wave state can be induced without 
the need for the external field. As previously established, 
s-wave order does not exist at half-filling as a mainfesta- 
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FIG. 4: Hysteresis of the superconducting order parameters. 
n s — |n s (K;)|. Starting from a d-wave state at half- 
filling, increasing the chemical potential increases the filling 
and decreases the d-wave order. Eventually, at n = 1.18 the 
system changes to an s-wave state. On return from large 
filling, the s-wave superconductivity is persistent to a low 
filling of n = 1.04, before spontaneously reverting to a d- 
wave state. The system is highly susceptible to d-wave order, 
and application of a very small external superconducting field 
to an s-wave state results in a d-wave state. Note that d- 
and s-wave channels are coupled in the higher order theory, 
so the transition can take place spontaneously, unlike in the 
standard gap equations. 



tion of Hohcnberg's theorem!, so the computed d-wave 
order at half-filling is the ground state of the model. It 
is interesting that the d- and s-channels are able to co- 
exist, considering that the BCS channels are separate on 
a square lattice. This is due to the vertex corrections, 
since the self consistent equations are no longer linear in 
the gap function (the 1st order gap equation vanishes in 
the d-wave case, leaving 2nd order terms as the leading 
contribution) . 

I finish with a brief discussion of Coulomb effects. 
In the Eliashberg equations, a Coulomb pseudopotential 
may be added to the theory as, 
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It is easy to see the effect of d-wave order on this term. 
Since the sign of the anomalous Green function is mod- 
ulated, the average effect of d-wave order is to nullify 
the Coulomb contribution to the anomalous self-energy 
(i.e. (j)cd = 0). This demonstrates that the d-wave 
state is stable to Coulomb perturbations, presumably be- 
cause the pairs are distance separated. In contrast, the 
s-wave state is not stable to Coulomb interaction, with 
a corresponding reduction of the transition temperature 
(Tc = for A < jUc)- Thus, such a Coulomb filter selects 
the d-wave state (see e.g. Ref. Ha) . Since large local 
Coulomb repulsions are present in the cuprates (and in- 
deed most transition metal oxides) , then this mechanism 
seems the most likely to remove the hysteresis. Without 
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the Coulomb interactions, it is expected that the s-wave 
state will dominate for n > 1.04, since the anomalous 
order is larger. 

I note that a further consequence of strong Coulomb 
repulsion is antiferromagnetism close to half-filling. Typ- 
ically magnetic fluctuations act to suppress phonon medi- 
ated superconducting order. As such, one might expect a 
suppression of superconducting order close to half-filling, 
with a maximum away from half filling. The current the- 
ory could be extended to include additional anomalous 
Green functions related to antifcrromagnetic order. This 
would lead to a 4x4 Green function matrix. A full anal- 
ysis of antiferromagnetism and the free energy will be 
carried out at a later date. 

a. Summary In this paper I have carried out simu- 
lations of the 2D Holstein model in the superconducting 
state. Vertex corrections and spatial fluctuations were 
included in the approximation for the self-energy. The 
anomalous self energy and superconducting order param- 
eter were calculated. Remarkably, stable superconduct- 
ing states with d-wave order were found at half-filling, 
d-wave states persist to n — 1.18, where the symmetry 
of the parameter changes to s-wave. Starting in the s- 
wave phase and reducing the filling, d-wave states spon- 
taneously appear at n = 1.04. The spontaneous appear- 
ance of (i-wave states in a model of electron-phonon in- 
teractions is of particular interest, since it may negate 
the need for novel pairing mechanisms in the cuprates 19 . 

The inclusion of vertex corrections and spatial fluctua- 
tions was essential to the emergence of the d-wave states 
in the Holstein model, which indicates why BCS and ME 
calculations do not predict this phenomenon. For very 
weak coupling, the off diagonal Eliashberg self-energy 
has the form — UT n F(iu) n , Q)Z?o(iw s — iuj n ), so it 
is clear (for the same reasons as the Coulomb pseudopo- 



tential) that this diagram has no contribution in the d- 
wave phase (the weak coupling phonon propagator is mo- 
mentum independent for the Holstein model). Therefore, 
vertex corrections are the leading term in the weak cou- 
pling limit. Furthermore, I have discussed the inclusion 
of Coulomb states to lowest order, which act to desta- 
bilize the s-wave states, while leaving the d-wave states 
unchanged. Since the Coulomb pseudopotential has no 
effect then it is possible that electron-phonon interactions 
are the mechanism inducing d-wave states in real mate- 
rials such as the cuprates. The Coulomb filtering mech- 
anism works for p-wave symmetry and higher, so it is 
possible that electron-phonon interactions could explain 
many novel superconductors. Certainly, such a mecha- 
nism cannot be ruled out. The doping dependence of the 
order qualitatively matches that of I^-aSr^CuO^ (here 
order extends to x = 0.18, in the Cuprate to x = 0.3). 
Antiferromagnetism is only present in the cuprate very 
close to half filling (up to approx x — 0.02), and on a 
mean-field level does not interfere with the d-wave su- 
perconductivity at larger dopings. 

It has been determined experimentally that strong 
electron-phonon interactions and high phonon frequen- 
cies are clearly visible in the electron and phonon band 
structures of the cuprates, and are therefore an essential 
part of the physics^. Similar effects to those observed 
in the cuprates are seen in the electron and phonon band 
structures of the 2D Holstein model in the normal phase 5 . 
It is clearly of interest to determine whether other fea- 
tures and effects in the cuprate superconductors could be 
explained with electron-phonon interactions alone. 
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